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ABSTRACT 

The large-scale field of the Sun is well represented by its lowest energy (or potential) 
state. Recent observations, by comparison, reveal that many solar-type stars show large-scale 
surface magnetic fields that are highly non-potential - that is, they have been stressed above 
their lowest-energy state. This non-potential component of the surface field is neglected by 
current stellar wind models. The aim of this paper is to determine its eff'ect on the coronal 
structure and wind. We use Zeeman-Doppler surface magnetograms of two stars - one with 
an almost potential, one with a non-potential surface field - to extrapolate a static model of 
the coronal structure for each star. We find that the stresses are carried almost exclusively in 
a band of uni-directional azimuthal field that is confined to mid-latitudes. Using this static 
solution as an initial state for an MHD wind model, we then find that the final state is de- 
termined primarily by the potential component of the surface magnetic field. The band of 
azimuthal field must be confined close to the stellar surface as it is not compatible with a 
steady-state wind. By artificially increasing the stellar rotation rate we demonstrate that the 
observed azimuthal fields can not be produced by the action of the wind but must be due to 
processes at or below the stellar surface. We conclude that the background winds of solar-like 
stars are largely unaff'ected by these highly-stressed surface fields. Nonetheless, the increased 
flare activity and associated coronal mass ejections that may be expected to accompany such 
highly-stressed fields may have a significant impact on any surrounding planets. 
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1 INTRODUCTION 

The magnetic fields of solar-like stars are an important influence 
not only on the rotational evolution of the stars themselves, but also 
on the atmospheres and exospheres of any planets that might sur- 
round them. This magnetic field not only transfers torques between 
the protoplanetary disk and the young star, but it also governs the 
loss of angular momentum in a wind. Any orbiting planets are ex- 
posed to the erosive effects of this wind and also the coronal X-ray 
emission from the star ( [Khodachenko et al.|2007| ). 

Both of these effects are likely to weaken as the star ages 
and spins down, generating less magnetic flux and hence produc- 
ing a weaker wind and reduced X-ray emission ( [Gudel|200 4). Re- 
cent maps of the surface magnetic fields of stars with a range of 
masses and rotation rates, however, suggest that it is not only the 
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strength of the magnetic field that changes with rotation rate, but 
also i ts geometry dPonati et al.|2008|porin et al.|2008l [Petit etlT] 
|2008||Morin et al.|2010| |. In contrast to the Sun which shows spots 
in well-defined "active latitudes", solar mass stars that are still in 
the rapidly-rotating stage typically show very non- solar magnetic 
fields, with spots that extend over the whole surface, often result- 
ing in a dark polar cap ( |Strassmeier|2009| ). Mixed polarity magnetic 
flux is seen at all latitudes on these stars. 

Typically these rapidly-rotating stars have X-ray luminosities 
that are three orders of magnitude greater than that of the Sun, but 
the extent of the corona that produces this emission is currently 
unknown. X-ray spectra suggest that their coronae are dense and 
compact (Dupree et al. 1993, Schrijver et al. 1995, Brickhouse &J 
Dupree 1998, Maggio et al. 2000, Gudel et al. 2001 , Sanz-Forcada] 
iet al.|20 03 ). In contrast, the presence of multiple large cool promi- 
nences trapped in co-rotation at distances of several stellar radii 
suggest that their closed magnetic fields, if not their X-ray bright 
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Figure 1. Surface magnetic field maps of CE Boo derived from spectropolarimetric observations ( [Donati et al.|2008) . The single black line shows the zero-field 
contour that separates regions of opposite polarity. 
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Figure 2. Surface magnetic field maps of GJ 49 derived from spectropolarimetric observations jDonati et al.|2008) . The single black line shows the zero-field 
contour that separates regions of opposite polarity. 



coronae, must be very extended (Collier Cameron & Robinson' 
1 989a b, [ Collier Cameron & Woods(r i992, Jeffries 1993, Byrne 
jetal. 1996, Elbe 1998, Barnes et al. 2000, Donati et al. 2000). One 
possible explanation is that the prominences form not within the X- 
ray bright corona, but in the cusps of helmet streamers that extend 
out into the stellar wind ( |Jardine & van Ballegooijen|2005| ). These 
prominences typically form in a timescale of 1 day and some 1-10 
appear in the observable hemisphere at any time. Their ejection in 
the stellar equivalent of solar coronal mass ejections not only con- 
tributes to the angular momentum loss from the star ( [Aarnio et al.| 
|2011a|b| ) but it will also temporarily enhance the ram pressure of 
the stellar wind and hence the degree of compression of any plane- 
tary magneto spheres. 

The surface fields of these young stars show one other very 
non-solar feature and that is the presence of a strong (some- 
times dominant) non-potential component ( [Petit etal. [2008 1 ). Stellar 
winds can of course produce azimuthal fields as the escaping wind 
extracts angular momentum from the star via magnetic torques, 
but for slow rotators it is unlikely that this could generate such 
strong fields at the photospheric level. Several other mechanisms 
have been proposed to explain the surface azimuthal fields, includ- 
ing the underlying dynamo ("Donati & Collier Cameron|1997| l, and 
the efifect of dififerential rotation in the presence of a unipolar cap 
( [Pointer et al.|2002| ). 

For solar mass stars, the surface dififerential rotation is similar 
to that of the Sun, but for higher-mass stars the dififerential rotation 
can be extreme, with equator to pole lap times as short as 16 days 
( [Barnes et al.|2005|[Marsden et al. [20051 [2006) [Jefers et al.|201l| ). 
The eflTect that this enhanced shear might have on the coronal and 
wind dynamics and the possible rate of coronal mass ejections is 



unknown. For low mass stars the dififerential rotation is typically 
weak i Morin et al.|2008| ). The high flaring rate of these stars how- 
ever suggests that some dynamic process is stressing the coronal 
field - even although in many cases the large-scale field that is de- 
tected with Zeeman-Doppler imaging is close to its potential, or 
lowest-energy state. 

Some insight into these stellar fields can be gained by consid- 
ering the changes in the solar magnetic field over the Sun's mag- 
netic cycle. At minimum, the solar field is closest to an aligned 
dipole, with fast wind streams emerging from the open field re- 
gions at the pole and the slow streams emerging from above the 
low-latitude active regions. As the cycle progresses, more bipoles 
emerge, contributing to the azimuthal field. These are acted on by 
dififusion and dififerential rotation, and their transport towards the 
poles by the meridional flow eventually reverses the polar polar- 
ity. In addition, their net contribution to the azimuthal field causes 
the axis of the large-scale dipole to move down into the equatorial 
plane and eventually reverse. This growth of active regions (and 
associated coronal mass ejections) through the cycle is also accom- 
panied by the extension of the polar coronal holes down towards the 
equatorial plane. As a result, fast wind streams originate at a range 
of latitudes and may interact with the slow wind streams to produce 
"co-rotating interaction regions" in the solar wind. These shocks 
provide a local density enhancement that, combined with the in- 
creased number of coronal mass ejections, can modulate the cosmic 
ray flux at Earth (Wang et al. 2006). Recent models of the variation 
of the solar wind through its cycle ( [Pinto et al.|2011| show that the 
magnetic torques exerted on the Sun vary significantly through its 
cycle, giving two orders of magnitude variation in the spin-down 
time. 



Table 1. Stellar and magnetic parameters for CE Boo and GJ 49, taken from 
[Donati et al. (2008 1. The table lists sequentially the stellar name, spectral 
type, mass, radius, inclination of the rotation axis and the rotation period, 
and then the field properties: the reconstructed magnetic flux density, and 
the fractional energy in the poloidal (potential) field. 
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CEBoo 
GJ49 


M2.5 
M1.5 


0.43 
0.51 


0.48 
0.57 


45 
45 


14.7 
18.6 


103 

27 


0.95 
0.48 



By analogy with the Sun, the very active young stars that show 
predominantly non-axisymmetric and non-potential surface fields 
may have winds that show a mixture of fast and slow wind streams 
with coronal mass ejections emerging from a range of latitudes. In- 
deed, the fact that these stars typically show mixed-polarity flux 
at all latitudes may suggest that their winds (while showing some 
characteristics of the solar wind at maximum) are much more ex- 
treme than the solar wind. 

Most stellar wind models are, however, based on the solar 
analogy. The simplest early models, such as the traditional Weber- 
Davies model ('Weber & Davies|1 967), assumed a split monopole, 
but more recent work usually initiates MHD simulations from an 
initial state defined by a "potential field source surface" model 
( |Altschuler & Newkirk, j"nl[T969l [Schatten et aLl[T969l ). This ap- 
proach assumes that the field is potential (i.e. in its lowest energy 
state) and that at some height above the surface the field lines are 
opened up by the pressure of the hot coronal gas. This method uses 
only the radial field component at the surface, neglecting the az- 
imuthal and meridional components. Its advantage is that it is com- 
putationally cheap and it provides a unique solution for the mag- 
netic structure. A recent comparison of the global structure pre- 
dicted by both the potential field source surface method and the 
full MHD simulation suggests that the former captures the large- 
scale structure of the solar coronal field fairly reliably ( [Riley et aT] 
|2006l l. 

This approach would not however capture the non-potential 
nature of the magnetic fields observed at the surfaces of other 
stars. The purpose of this paper is to explore the eflfect of this non- 
potential field on the large-scale structure of the corona and winds 
of solar- type stars. 



2 THE SURFACE MAGNETOGRAMS 

In order to study the eflfects of the non-potential field on the coronal 
structure and dynamics we choose to compare two stars (CE Boo 
and GJ 49) that are similar in rotation rate but with slightly diflfer- 
ent masses. One has a surface field that is close to potential, while 
the other has a significant non-potential component. Both stars are 
slow rotators, so rotational eflfects are minimal. In addition, the in- 
clination of the rotation axes of both stars to the line of sight is 
the same, so the magnetic fields of both stars are seen in the same 
orientation. The stellar parameters are shown in Table [T] 

We choose initially to compare the static coronal structures 
that are found by assuming either that the field is purely potential, 
or that it has both potential and non-potential components. These 
extrapolations can be used as the initial condition for a full MHD 
solution. Since we are particularly interested in the non-potential 
field, we also explore the possibility that for the more rapidly- 



rotating stars it is the rotational stressing of the surface field by 
the action of the wind that causes the field to depart from a poten- 
tial state. We therefore perform one simulation of the wind of GJ 
49 with an artificially decreased stellar rotation period of 0.6 days 
and compare this to the wind parameters found with the observed 
rotation period of 18.6 days. 

The input for the static extrapolation is taken from |Donatietal.| 
([2008). The surface magnetic field of both stars were modelled with 
Zeeman-Doppler imaging from time- series of spectropolarimetric 
observations collected over approximately 2 consecutive stellar ro- 
tations. For spatially unresolved sources, due to the mutual cancel- 
lation of contributions from neighbouring regions of opposite po- 
larities to the polarized signal, spectropolarimetric measurements 
can only probe the large-scale component of magnetic fields (see 
e.g. [Morin| ( [2012| )). The maximum degree i of modes that can be 
reconstructed with ZDI depends on the star's projected rotational 
velocity. For slow rotators such as CE Boo and GJ 49, the recon- 
struction is limited to modes with order i <, 8. As there is no 
unique solution to the ZDI problem, a regularization scheme has 
to be used. A maximum entropy solution corresponding to the low- 
est magnetic energy content is used. It is optimal in the sense that 
any feature present in the map is actually required to fit the data. 
Although this method does not allow us to derive formal error bars 
on the reconstructed maps, numerical experiments have shown that 
ZDI is a robust method ( [Donati & Brown|1997[[Morin et al.|20"T0l ). 

This reconstructed field is expressed as a sum of a poloidal 
and toroidal field ( [Mestel|1999| ). The poloidal component captures 
the potential contribution to the total field, that is the component 
that is in its lowest-energy state. The toroidal component lies on 
the surfaces of concentric spheres and captures the non-potential 
component of the total field. It is this component that is associ- 
ated with the electric currents in the corona and which describe the 
free energy that is available to power, for example, stellar flares 
and coronal mass ejections. These two components of the surface 
field can be expressed as linear combinations of spherical harmon- 
ics ( [Donati et~al. 2006 ). Thus the radial, meridional and azimuthal 
field components at the stellar surface can be written in spherical 
co-ordinates (r, 6, 0) aQ 
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and Pim(0) denotes the associated Legendre functions. The poten- 
tial terms are those with coefficients ai,n or fii^, while the non- 
potential terms are those with coefficients yi^. Clearly, then, in the 
limit yim ^ we recover a purely potential field. 



^ We note that in [Donati et al.[{2006) , the radial field is positive outwards, 
the azimuthal field is positive in the direction of stellar rotation (i.e. in- 
creasing longitude or decreasing rotation phase) and the meridional field is 
positive when pointing to the visible pole. 
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Figure 3. The static solution for the surface magnetic field of CE Boo, divided into its diff'erent components. The meridional component is shown in the left 
column and the azimuthal component in the right column. The top row shows the non-potential contribution and the bottom row the potential contribution to 
the total field. The single black line shows the zero-field contour which therefore separates regions of opposite polarity. 
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Figure 4. The wind solution for the surface magnetic field of CE Boo, divided into its meridional (left column) and azimuthal (right column) components. The 
single black line shows the zero-field contour which therefore separates regions of opposite polarity. 



The corresponding surface magnetic maps from |Donati et al.| 
( |2008| ) are reproduced in Figs. [T] and [2] In both cases, the radial 
and meridional fields look very similar to a dipole, but particularly 
in the case of GJ 49, there is a significant azimuthal field that is 
unidirectional at low to mid latitudes. This is a clear signature of a 
non-potential field. 



3 THE STATIC CORONAL MAGNETIC FIELD 

In order to determine the coronal structure that corresponds to these 
surface fields we need to make some assumptions about the nature 
of the coronal field. The simplest assumption is that the field is po- 
tential, or in its lowest energy state and is determined simply by the 
coefficients aim andy^/^ in Eqns. ([l]-[3]). This is the starting point 
for many extrapolations of the solar magnetic field. If we wish to 
determine the distribution of electric currents in the corona, how- 
ever, we need to allow for the non-potential components that are 
described by the coefficients y/^. 



3.1 Potential field extrapolation 

We begin by calculating the contribution to the total field that is 
potential. We write B^^^ in terms of a ffux function ^ such that 
^pot ^ _ yxj, ^^^ ^j^g condition that the field is potential (V X 5^"^^ = 
0) is then satisfied automatically. The condition that the field is 
divergence-free then reduces to Laplace's equation V^^ = with 
solution in spherical co-ordinates (r, 6, 0) 



^-J]J][aMr^ + bimr-^^^'^ 



]Pim{e)e'' 



(5) 



where all radii are scaled to a stellar radius and the associated Leg- 
endre functions are once again denoted by P/^. The two unknowns 
are therefore the coefficients aim and bim- One of these can be deter- 
mined by imposing the radial field at the surface from the Zeeman- 
Doppler maps. In order to determine the second unknown, we select 
a particular form of potential field that has the useful property that 
at some radius, all the field lines are open. This mimics the eff'ect of 
the outward pressure of the hot coronal gas pulling open field lines 
to form the stellar wind. Thus, at some normalised radius Rs above 
the surface (known as the source surface) the field becomes radial 
and hence Bo(Rs) = ^^(P J = 0. As a result. 
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where the functions /;(r, R,) and gi(r, Rs) which describe the influ- 
ence of the source surface (and hence the wind) on the magnetic 
field structure are given by 



Mr,Rs) = 
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/+l+/(l/i?,)2'+i 



(10) 



gi(r,Rs)- 



1 - (r/Rsf 



l+/(l/i?,)2'-n 



(11) 



In the limit where the source surface is large (i.e. the magnetic field 
is completely closed), we recover the familiar multipolar expan- 
sions for a magnetic field. This limit corresponds io Rs ^ oo and 
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The coefficients Bim are determined by the surface radial field that 
is derived from the Zeeman-Doppler maps (i.e. by the values of a^ 
in ITJ). This is known as the Potential Field Source Surface method. 
It was originally developed for extrapolating the Sun's coronal field 
from solar magnetograms ( [Altsc huler & Newkirk, Jr. 1969 ). We use 
a code originally developed by van Ballegooijen et al., ( ,1998| ) (see 
also |Jardine etaT] ( [20021 )). 

Comparing the form of our extrapolated field given in ([7] -[9]) 
with the general expressions for the observed field at the surface 
([T]-[3]) we can see that our extrapolated field cannot match the ob- 
served surface field exactly. The reason is that the form of potential 
field we are using for the extrapolation (the Potential Field Source 
Surface method) is only one type of potential field. The assumption 
of a source surface forces a relationship between the field compo- 
nents that means they are no longer independent. While aim can be 
simply related to Bim, we cannot match the values of /?/,„ that are de- 
rived from the observations. Therefore this method, which selects 
only one type of potential field, will not be guaranteed to reproduce 
the potential field contribution to Bq and B^ that is fitted to the data. 

With this caveat in mind, we use the observed By at the stellar 
surface to determine Bim and hence to obtain the potential contribu- 
tion to the azimuthal and meridional fields B^^^ and B^^^ We show 
these in the bottom rows of Fig.[3]and Fig.|5] It is clear by compar- 
ison with the observed surface maps shown in Figs. [T] and |2] that 
this potential field does not reproduce all the observed field com- 
ponents. In particular, the uni-directional band of azimuthal field is 
absent from these potential field maps. In order to extrapolate the 
non-potential part of the field, however, we need to make an as- 
sumption about the nature of the coronal currents. We base our ex- 
trapolation on the method developed by |Hussain et al.|p002| ). This 
is not a force-free solution, but it allows us to incorporate fully the 
non-potential contribution of the surface field and to extrapolate it 
into the corona. 



3.2 Non-potential field extrapolation 

In general, the magnetic field will be a sum of potential and non- 
potential terms such that 5 = 5^^^ + 5"^. We assume that the non- 
potential magnetic field is perpendicular to the radial direction (i.e. 
it lies on spherical shells and so B^/ - 0). Furthermore, the electric 
currents are assumed to be derived from a potential Q: 



V X B"P = -VQ. 



(14) 



It follows that V^2 = 0, so g(r) has a solution in terms of spherical 
harmonics. As shown in the Appendix, we find solutions for this 
non-potential magnetic field that vanish at the source surface and 
have the form 



5;^ = o 



(15) 
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Figure 5. The static solution for the surface magnetic field of GJ 49, divided into its diff'erent components. The meridional component is shown in the left 
column and the azimuthal component in the right column. The top row shows the non-potential contribution and the bottom row the potential contribution to 
the total field. The single black line shows the zero-field contour which therefore separates regions of opposite polarity. 
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and disRs ^ oo we recover /z/(l) -^ \/l. 

While this is not the most general form of non-potential field, 
it has the useful property that the equations for B^^^ and 5"^ are 
now structurally very similar to the forms used in ([l]-|3]) to describe 
the surface field. The coefficients yim and Cim that govern the non- 
potential field components can be simply related. This allows us 
to match the observed non-potential component of the field exactly 
to our model and to extrapolate it into the corona. Thus while the 
potential part of our extrapolated field will not reproduce an exact 
match to the potential part of the observed surface field, the non- 
potential part matches exactly. 

We therefore show the non-potential (top row) and potential 
(bottom row) parts of the field separately in Figs. [3] and |5] The to- 
tal field is the sum of both of these. Fig. |7] shows the extrapolation 
of this total field with a source surface chosen to be at 4R^. The 
largest closed field lines have been selected in order to highlight 
the structure of the large-scale field. The tilt of the dipole axis can 
be clearly seen in both cases, although it should be noted that the ro- 
tation axes of both stars have the same inclination to the observer's 
line of sight. While the extrapolation of the potential contribution to 



the total field is fairly similar in both stars, the inclusion of the non- 
potential contribution highlights the dififerences between the mag- 
netic field structures of the two stars. The non-potential component 
introduces an azimuthal shear into the field that is most apparent in 
GJ 49 (for which 52% of the total magnetic energy in the surface 
field is non-potential). 



4 THE STELLAR WIND MODEL 

To perform the stellar wind simulations, we use the three- 
dimensional MHD numerical code BATS-R-US developed at Uni- 
versity of Michigan (Pow ell et al.| „1999). BATS-R-US has been 
widely used to simulate, e.g., the Earth's magnetosphere (Rid-] 
[ley et al. 2006), the heliosphere (R oussev et a l. 2003), the outer- 
heliosphere (Linde et al.'r T998nOpher et aL||2003,|2004) , coronal 
mass ejections (Manchester et al. |2004| |Lugaz et al.||2b 05 ), the 
magnetosphere of planets ( Toth et al. 2004 , Hansen et al.||2005] l, 
and stellar winds of cool stars (Vidotto et al. 2009 2012). It solves 



the ideal MHD equations, that in the conservative form are given 
by 



dp 



^ + V • (pu) = 0, 




dt 




BB 

47T 



= PS, 



(19) 
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Figure 6. The wind solution for the surface magnetic field of GJ49, divided into its meridional (left column) and azimuthal (right column) components. The 
top row shows the result of assuming the observed stellar rotation period of 18.6 days, while the bottom row shows the result of assuming a stellar rotation 
period artificially decreased to 0.6 days. The single black line shows the zero-field contour which therefore separates regions of opposite polarity. 



'dt 
ds 

'di 



+ V (uB-Bu) = 0, 

B2\ (u B)B 



+ V 



u e + p + 



8;r 



4;r 



= Pg u, 



(21) 

(22) 



where the eight primary variables are the mass density p, the plasma 
velocity u = [Ur, Ug, u^], the magnetic field B = [Br, Bq, 5^}, and the 
gas pressure p. The gravitational acceleration due to the star with 
mass M^ and radius R^ is given by g, and s is the total energy 
density given by 



pu p 

2 7- 



1 "^ 8;r' 



(23) 



We consider an ideal gas, so p = nksT, where ks is the Boltzmann 
constant, T is the temperature, n = p/iprrip) is the particle number 
density of the stellar wind, finip is the mean mass of the particle, 
and 7 is the poly tropic index (such that p oc p^). 

As the initial state of the simulations, we assume that the wind 
is thermally driven { [Parker|1958| ). At the base of the corona (r = 
R^),WQ adopt a wind coronal temperature Tq = 2 x 10^ K and wind 
number density hq = lO^^cm"^. The stellar rotation period Prou 
M^ and R^, are given in Table [T] With this numerical setting, the 
initial solution for the density, pressure (or temperature) and wind 
velocity profiles are fully specified. 

To complete our initial numerical set up, we assume that the 
magnetic field is either potential everywhere (i.e., V X B = 0) or the 
sum of potential plus non-potential components, as described in 



sections ( |3.1| ) and ( |3.2| ). The initial solution for B is found once the 
distance to the source surface is assumed (set at 4 R^ in the initial 
state of our runs) and the surface magnetic field is specified: either 
simply the radial component (in the case of a potential field) or all 
three components (in the case of a total potential plus non-potential 
field). 

Once set at the initial state of the simulation, the distribution of 
Br is held fixed at the surface of the star throughout the simulation 
run, as are the coronal base density and thermal pressure. A zero 
radial gradient is set to the remaining components of B and u = 
in the frame corotating with the star. The outer boundaries at the 
edges of the grid have outflow conditions, i.e., a zero gradient is set 
to all the primary variables. The rotation axis of the star is aligned 
with the z-axis, and the star is assumed to rotate as a solid body. 

Our grid is Cartesian and extends in x, y, and z from -20 to 
20 R^, with the star placed at the origin of the grid. BATS-R-US 
uses block adaptive mesh refinement, which allows for variation in 
numerical resolution within the computational domain. The finest 
resolved cells are located close to the star (for r <2R^), where the 
linear size of the cubic cell is 0.02 R^. The coarsest cell is about 
one order of magnitude larger (linear size of 0.31 R^) and is lo- 
cated at the outer edges of the grid. The total number of cells in our 
simulations is about 15 million. 

As the simulations evolve in time, both the wind and magnetic 
field lines are allowed to interact with each other. The resultant so- 





Figure 7. Static field line extrapolations for CE Boo (top) and GJ 49 (bottom) for fields that are purely potential (left) and those that are the sum of potential 
plus non-potential (right). Closed field lines which would contain coronal gas are shown white, open field lines which would contribute to the stellar wind are 
shown blue. 



lution, obtained self-consistently, is found when the system reaches 
a steady state (in the reference frame corotating with the star). 



5 RESULTS AND DISCUSSION 

We have separated the magnetic fields of CE Boo and GJ 49 into 
their lowest-energy (potential) and stressed (non-potential) com- 
ponents. This has allowed us to isolate both the locations where 
the field is stressed above its lowest energy state and also the na- 
ture of the structures that carry these stresses. We find that the de- 
partures from a lowest-energy state are apparent mainly in the az- 
imuthal field (the meridional field contributes a negligibly- small 



non-potential component). This appears as a clearly-defined mid- 
latitude band of unidirectional azimuthal field (see Figs. [3] and [Sj. 
This is similar to the non-potential field of the young rapid rota- 
tor AB Dor { [Hussain et al. 2002) except that it appears at lower 
latitudes. 

By extrapolating these surface fields into the corona we can 
see that the presence of the non-potential field does not change the 
overall topology of the coronal field, but it provides an azimuthal 
shear (see Fig. [7]). We have also explored the nature of the winds 
that might be associated with these surface fields by using these ex- 
trapolated fields as an initial state for an MHD wind model. As the 
solution evolves towards a steady- state, only the radial component 



CEBoo 




of the surface field is kept fixed - the azimuthal and meridional field 
components are allowed to vary in response to the available forces. 

The top and middle rows of Fig. IS] show the final structure of 
the magnetic fields of CE Boo and GJ 49. This final state was the 
same, regardless of whether the initial state was the total field (ie the 
potential plus non-potential field), or simply the potential field. This 
happens because it is the boundary conditions that will control the 
final state and the initial conditions of the system will be "flushed 
out" by the wind. Therefore, we find that the mass loss rates, angu- 
lar momentum loss rates, and the fluxes of surface magnetic field 
and open magnetic field are the same in both cases. Since there are 
no forces in the wind model capable of providing the stresses nec- 
essary to sustain the strong azimuthal field at the stellar surface, the 
solution relaxes back to something close to a potential field at the 
surface. As a result, the steady-state wind solution has meridional 
and azimuthal field components at the stellar surface that are close 
to the potential component of the field that we calculate and that do 
not reproduce the observed non-potential component of the surface 
field. 

At larger heights of course, the action of the wind stresses the 
field and generates an azimuthal field component, but this is fairly 
small at the stellar surface, particularly for such slowly-rotating 
stars. In order to confirm the role of rotation in influencing the field 
structure, we also artificially increase the rotation rate of GJ 49, 
while keeping the initial magnetic field structure unchanged. The 
resulting field structure is shown in Fig.lS] While an azimuthal field 
develops with height in the corona, it is small at the surface and can- 
not explain the observations. This suggests that this azimuthal field 
is produced not by the wind, but by the sub-surface dynamo. 

These simulations therefore suggest that the ambient winds of 
these slowly-rotating stars are well described by the potential com- 
ponents of their surface fields. The strong azimuthal fields seen at 
the surface should not survive to the heights in the corona at which 
the wind is launched. They may of course be important in deter- 
mining flare locations and energies. For GJ 49, for example, 52% 
of the total magnetic energy close to the surface is contained in the 
non-potential part of the field and is therefore available for release. 
It is mainly contained in a well-defined band that is centred around 
latitudes 30° - 40°. This is the maximum latitude at which solar ac- 
tive regions are seen and from which solar coronal mass ejections 
are launched. This might suggest that this is the region from which 
flares and coronal mass ejections could be expected. On the young 
rapid rotator AB Dor, by comparison, (Prot = 0.514 days) the band 
of non-potential field is strongest around 70° - 80° ( [Hussain et al.| 
|2002| ) which may suggest a diff'erent pattern of coronal mass ejec- 
tion. Such coronal mass ejections would temporarily increase the 
mass loading of the stellar wind and also its ram pressure, which 
is responsible for compressing the magnetospheres of any orbit- 
ing planets. Whether coronal mass ejections provide a significant 
contribution to either angular momentum loss or the impact of the 
wind on orbiting planets depends on their size and frequency. The 
background stellar wind that we find however is independent of the 
strong non-potential component of the surface fields and is primar- 
ily governed by their radial component. 



Figure 8. Final magnetic field structures for CE Boo (top row) and GJ 49 
(middle row). The final state is the same, regardless of whether the initial 
state is the total field (i.e. the potential plus non-potential field), or simply 
the potential field. The bottom row shows the eff'ect on GJ 49 of artificially 
decreasing the rotation period from 18.6 to 0.6 days 



APPENDIX A: THE NON-POTENTIAL FIELD 

We look for solutions for the non-potential field that are of the form 
\ dF \dF 



Bm 



■ 0, 5n, 



r sin ^ < 



B, 



np,0 ■ 



r do' 



(Al) 
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where F(r) is a scalar function. These automatically satisfy V • 5 = 
0. Furthermore, the electric currents are assumed to be derived from 
a potential Q: 



VxBn 



-vg. 



(A2) 



It follows that V^2 = 0, so Q(r) has a solution in terms of spheri- 
cal harmonics. Inserting equation \A\\ into ( |A2| ), we find from the 
radial component: 



i^l""^ 



do 



1 ^2^ 






and from the ^ and components: 



'rdr\~de 



IdQ I d I I dF 

r 86 r dr\ sin 6 d(p 



_\_dQ 
r sin 6 dd) 



which can be integrated with respect to 6 and 0: 
dp 

We now introduce a third scalar C(r) such that 

"^' 
then equation (|A5|) yields 



F = r^ 



0--—1 ^ — 
dr \ dr 



(A3) 



, (A4) 



(A5) 



(A6) 



(A7) 



Inserting these expressions for F and Q into equation ( |A3| ) we ob- 



tam: 



(r^V^c) = 0. 



(A8) 



Assuming V^C = at the stellar surface, it follows that this condi- 
tion is true at all heights, so C(r) is also a harmonic function. We 
now write C as a sum over spherical harmonics: 



C{r, e,^) = )^ C,„q,{r)P,,Me""*, 

Im 



(A9) 



where / is the harmonic degree (/ = 1,2, • • •), w is the azimuthal 
mode number (-1 ^ m ^ +/), Pim(0) is the associate Legendre 
function, and qi(r) describes the radial dependence of the various 
modes. Note that the function qi(r) defined here is diff'erent from 
the function //(r) for the potential field (see Hussain et al. 2002). 
Then the non-potential components of the magnetic field follow 
from equation ( |A6| ) and \A\\ : 



^ dr smO 

Im 



(AlO) 



(All) 



The function qi(r) must satisfy the following constraints. First, we 
assume that qi(l) = 1, so that C^ are the mode amplitudes at the 
stellar surface (r = 1). Second, since C(r) is a harmonic function, 
qi(r) must be a sum of a radially decreasing term oc r~^~^ and an 
increasing term oc rK Third, we require that the horizontal compo- 
nents of the non-potential field vanish at the source surface; this 
implies dqi/dr = at r = /?s. From these conditions it follows that: 



,^=-/(/+l)r<'-%(r,i?,) 
dr 



(A12) 



and hence we recover expressions ( 16 1 and ( 17 1 for Bf and Bj 
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